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The strategy and techniques for analyzing anisotropic flow (directed, 

ON \ 

ON . elliptic, etc.) in relativistic nuclear collisions are presented. The emphasis is 

£^ , on the use of the Fourier expansion of azimuthal distributions. We present 

formulae relevant for this approach, and in particular, show how the event 
' multiplicity enters into the event plane resolution. We also discuss the role of 

non-flow correlations and a method for introducing flow into a simulation. 
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I. INTRODUCTION: SUBJECT AND TERMINOLOGY 



^ , Recently, the study of collective flow in nuclear collisions at high energies has attracted 

increased attention of both theoreticians and experimentalists |flj . There are several reasons 



for this: i) the observation of anisotropic flow at the AGS H-H] an d at the SPS PHTT 
ii) progress in the theoretical understanding of the relation between the appearance and 
development of flow during the collision evolution, and processes such as thermalization, 
creation of the quark-gluon plasma, phase transitions, etc. [P^TSj? iii) the study of mean 



field effects [14,19], iv) the importance of flow for other measurements such as identical 



and non-identical two-particle correlation analyses |20|-42"3"|| , v) the development of new tech- 
niques suitable for flow studies at high energies P3-EBl|I],P[ . Although all forms of flow are 
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interrelated and represent only different parts of one global picture, usually people discuss 
different forms of collective flow, such as longitudinal expansion, radial transverse flow, and 
anisotropic transverse flow of which the most well established are directed flow and ellip- 



tic f2|,[19l flow [0. 

The study of flow in ultrarelativistic nuclear collisions has just begun and it is very 
important to define the subject and establish the terminology. At high energies the longi- 
tudinal flow is well decoupled from transverse flow. The flow (polar) angles observed at low 
energies are relatively large and a rotation of the coordinate system was done in order to 
analyze the event shape in the plane perpendicular to the main axis of the flow ellipsoid. 
At high energies the flow angles become very small, 9fi ow ~ (p x )/Pbeam <S 1, so that not 
doing such a rotation "induces" an elliptic anisotropy of the order of the square of the flow 
angle and usually can be neglected. This means that at high energies one does not have 



to rotate to the flow axis to study the flow pattern [27|, but one can use the plane trans- 



verse to the beam axis. Thus we discuss only anisotropic transverse flow from the particle 
azimuthal distributions at fixed rapidity or pseudorapidity. It appears very convenient to 



describe the azimuthal distributions by means of a Fourier expansion ||2q , |2q , |24| , and below 
we characterize different kinds of anisotropies as corresponding to different harmonics (in 
analogy/contrast to the description at low energies, where the three-dimensional event shape 
is characterized using multipole terminology). Anisotropic flow corresponding to the first 
two harmonics plays a very important role and we use special terms for them: directed and 
elliptic flow, respectively. The word "directed" comes from the fact that such flow has a 
direction, and the word "elliptic" is due to the fact that in polar coordinates the azimuthal 
distribution with non zero second harmonic represents an ellipse. 

II. CORRELATIONS WITH RESPECT TO AN EVENT PLANE 

Strategy. In this section we summarize an approach to the study of anisotropic flow 
which is particularly suitable at high (AGS/SPS/RHIC) energies. It uses the Fourier ex- 
pansion of azimuthal distributions, introduced in this way for the analysis in ref. f2"5fl . The 



essence of the method is to first estimate the reaction plane. The estimated reaction plane 
we call the event plane. The Fourier coefficients in the expansion of the azimuthal distribu- 
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tion of particles with respect to this event plane are evaluated. Because the finite number 
of detected particles produces limited resolution in the angle of the measured event plane, 
these coefficients must be corrected up to what they would be relative to the real reaction 
plane. This is done by dividing the observed coefficients by the event plane resolution, which 
is estimated from the correlation of the planes of independent sub-events (sub-groups 
of the particles used for the event plane determination). The resolution obtained from the 
sub-events can be converted to that for the full event by means of the multiplicity depen- 
dence of the resolution which will be described below. Also, if the detector does not have 
full azimuthal acceptance, the acceptance bias has to be removed. 

Fourier expansion. The quantity under study in the most general case is the triple 
differential distribution. In this, the dependence on the particle emission azimuthal angle 
measured with respect to the reaction plane can be written in a form of Fourier series 

where \l/ r denotes the (true) reaction plane angle, and the sine terms vanish due to the reflec- 
tion symmetry with respect to the reaction plane. The main advantage of the Fourier method 
is that the Fourier coefficients, evaluated using observed event planes, can be corrected for 
the event plane resolution caused by the finite multiplicity of the events. This correction 
always raises the value of the coefficients. The great importance of this is that then the 
results for particles in a certain phase space region may be compared directly to theoretical 
predictions, or to simulations unfiltered for the detector acceptance, and for which the reac- 
tion plane has been taken to be the plane containing the theoretical impact parameter. Note 
the factor of "2" in front of each v n coefficient. We propose to use it because in this case 
the meaning of the coefficients v n becomes transparent P5 |, v n = (cos(n(0 — ^ r ))), where 



( ) indicates an average over all particles in all events. For the particle number distribution, 
the coefficient v x is (p x /p t ) and v 2 is ({p x /pt) 2 - iPy/Pt) 2 )- 

Estimation of the reaction plane. The method uses the anisotropic flow itself to 
determine the event plane. It also means that the event plane can be determined indepen- 
dently for each harmonic of the anisotropic flow. The event flow vector Q n and the event 
plane angle from the n— th harmonic of the distribution are defined by the equations 
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Q n cos(n^ n ) = X n = w i cos(n0i) , (2) 

i 

Q n sm(n^ n ) = Y n = J2w i sm(n(j) i ), (3) 



or 



*n = tan — — — n. 4 

V Ei Wi cos(n0i) J 

The sums go over the i particles used in the event plane determination and the Wi are 
weights. In general the weights are also optimized to make the reaction plane resolution the 
best that is possible. Sometimes it can be done by selecting the particles of one particular 
type, or weighting with transverse momentum of the particles, etc. Usually the weights for 
the odd and even harmonic planes are different. Optimal weights are discussed in footnote 2 
of ||17|| . For symmetric collisions reflection symmetry says that particle distributions in the 
backward hemisphere of the center of mass should be the same as in the forward hemisphere 
if the azimuthal angles of all particles are shifted by 7r. This explains why for the odd 
harmonics the signs of the weights are reversed in the backward hemisphere while for the 
even harmonics the signs of the weights are not reversed. Note that the event plane angle 
\l/ n determined from the n— th harmonic is in the range < < 2ir/n. For the case of 
n — 1, Eqs. (|2||||4]) are equivalent to obtaining ^ for number flow from [2£ 



Q = J2 w Pt/\Pt\ (5) 
where the sum is over all the particles. The case of n = 2 is equivalent to the event plane 



determined from the transverse sphericity matrix |24 . 

Acceptance correlations. Biases due to the finite acceptance of the detector which 
cause the particles to be azimuthally anisotropic in the laboratory system can be removed 
by making the distribution of event planes isotropic in the laboratory. We know of several 
different methods to remove the effects of anisotropy which have been used (sometimes in a 
combination with each other). Each of them has some advantages along with disadvantages. 

The simplest one is to re-center P,|^ JT0|j30[| the distributions (X n , Y n ) (Eqs. |||3]) by sub- 



tracting the (X n , Y n ) values averaged over all events. The main disadvantage of this method 
is that it does not remove higher harmonics from the resulting distribution of \l/ ra . If such 
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harmonics are present then the method requires additional flattening of the event plane dis- 
tribution by one of the other methods. The second, one of the most commonly used methods, 
is to use the distribution of the particles themselves as a measure of the acceptance HdH. 
One accumulates the laboratory azimuthal distribution of the particles for all events and 
uses the inverse of this as weights in the above calculation of the event planes. The limitation 
of this approach is that it does not take into account the multiplicity fluctuations around 
the mean value. The third method is to use mixed events [|3],^||| . Correlations with the 
raw event planes are stored in histograms and correlations with event planes from previous 
events are also stored. The real correlations are then divided by the mixed event correlations 
to remove the acceptance correlations. This third method suffers from the problem that, 
if one uses only one mixed event for each real event, the errors are v2 larger. If one uses 
many mixed events for each real event, the errors decrease as n l J^ x instead of ^/n mix because 



the same events are being used n mix times [31]. The fourth method fits the unweighted 



laboratory distribution of the event planes, summed over all events, to a Fourier expansion 
and devises an event-by-event shifting of the planes needed to make the final distribution 
isotropic 0,|J[]- In all these methods one has to check that the event plane distributions 
indeed become isotropic. 

Particle distributions with respect to the event plane. The next step is to study 
the particle distributions with respect to the event planes. Note that for a given n the 
corresponding Fourier coefficient v n can be evaluated using the reaction planes determined 
from any harmonic m, with n > m, if n is a multiple of m. If n > m, the sign of v n is 
determined relative to v m . That is, the first harmonic plane can be used, in principle, to 
evaluate all v n . The second harmonic plane can be used to evaluate t> 2 , ^4> etc. For the 
event plane evaluated from the m-th harmonic the Fourier expansion is 



1 The equation for the shift is (see Appendix in ||] 

'&max ^ 

We have usually taken i max = 4/n for n = 1, 2 



A^ n = ^2 -(-(sin(m*„)) cos(m* n ) + (cos(in* n )) sin(in^ n )). (6) 
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d{wN) {wN} H + f>£>s(M0 -*-)))■ (7) 



Writing the equation in terms of km instead of n insures, for instance, that when m — 1 
all terms are present, but when m = 2 only the even terms are present. The quantity w 
is a weight which could be p t of the particle if one studies transverse momentum flow, or 
just unity, if one studies flow of the number of particles. The coefficients v are evaluated 
by (cos(n(0 — ty m ))). The quantity (0 — \l/ m ) has a lowest order periodicity of 2n/m. To 
graph the distribution one can shift the negative values to positive by adding 2ir, and then 
fold the distribution with this periodicity using the module function. When a particle has 
been used in the calculation of an event plane, the auto-correlation effect in its distribution 
with respect to this plane is removed by recalculating that plane without this particle []29 



This is easily done if one saves the sums of sines and cosines from Eqs. (0-^j), and subtracts 
the contribution of the particle and its weight from these sums. This method of removing 
autocorrelations assumes that contributions from conservation of momentum are small. 

The event plane resolution. The coefficients in the Fourier expansion of the azimuthal 
distributions with respect to the real reaction plane are then evaluated by dividing by the 
event plane resolution ^|^J2||||,[§], 



v n = v° n bs /(cos(km(y m -y r )). (8) 

The mean cosine values are less than one and thus this correction always increases the flow 
coefficients. 

The resolution depends both on the harmonic used to determine the event plane m and 
the order of the calculated coefficient n. It is generally true that better accuracy for the 
determination of v n is achieved by using the event plane (^ n ) determined from the same 
harmonic {m = n, k = 1) because the resolution deteriorates as k increases (see below). 
For example, better accuracy for v 2 can be achieved using \l/ 2 even when the elliptic flow is 
somewhat weaker than the directed flow. 

To calculate the resolution we start with the distribution of m(^ m — ^ r ), which can be 



written as |?5] 



dP _ f v' m dv' m . v^ + vZ- ^mv'm cos(m(^ m - fr r )) 



d(m(tf m -tf r )) J 2na 2 ^ v 2a 2 



The parameter a, which to second order in flow is common for all m, is inversely proportional 
to the square-root of N, the number of particles used to determine the event plane 



a 



2 



1 (w 2 ) 



(10) 



2N (w) 2 ' 

The integral @ can be evaluated analytically [p4,25[] , and then the event plane resolution 
can be expressed as 

(cos(fcro(*m - * P ))) = ^= Xm exp(-x^/4) [I^^Ji) + J^i(x^/4)] , (11) 

where Xm = v m /a (= v m \ / 2N for number flow)0 and I v is the modified Bessel function of 
order v. This resolution function is plotted in Fig. 1. 

Note that (cos(m(\E' jn — \I/ r ))) is a correction for the reconstruction of the signal if the 
event plane is derived using the flow of the same harmonic (k = 1). In the case when 
the harmonic orders do not coincide, for example, when one uses the event plane derived 
from the first harmonic (directed flow) and studies the second harmonic (elliptic flow) of 
the particle distribution with respect to this plane (the m = 1, k = 2 term in Eq. (^)) the 
correction would be (cos(2(\E' 1 — \P r ))). For practical use all such functions can be calculated 
numerically as in the section on Approximations below. 



The resolution correction used in earlier times at the Bevalac [28] is close to the present 
curve for values of xi = V\y2N greater than about 1.5, which was generally true at the 
Bevalac. However, using the old procedure at the higher beam energies of the AGS or SPS 
would greatly overestimate the resolution and make the flow values too small. 

The correlation between flow angles of independent sets of particles. If one 
constructs the event planes in two different windows, (a) and (b), or from two random sub- 
events0, the corresponding correlation function also can be written analytically. But more 
important in this case is the simple relation for such correlations, 



2 Please note that the parameter x used in this paper is a factor of v2 larger than the one used 
in refs. Iffi^]; it is equivalent to the parameter £ in [25|. 

3 The random sub-events can be made to have equal multiplicity using the CERN library routine 
SORTZV. 
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(cos(n(¥£ - *i))> = (cos(n(^ - ¥ r )))<cos(n(*i - ¥ P ))). ( 12 ) 

Here, the assumption is made that there are no other correlations except the ones due to flow, 
or that such other correlations can be neglected. (For example, two-particle correlations due 
to resonance decays should scale with multiplicity as 1/N and usually can be neglected.) 
If this is not true, then special precautions have to be made to avoid or correct for such 
correlations. See the section below on Non-flow Correlations. 

Note that the correlation of two planes (the distribution in \l> a — \|> fe ) is not well represented 
by a Fourier expansion. For a strong correlation it should be Gaussian. For no flow the 
correlation between two random planes of the same order should be a triangular distribution. 
For plotting purposes, if one takes the absolute value of — ^ b n and folds this distribution 
about the angle n/n, the triangle becomes flat. This does not affect the (cos(n( 1 i r * — 
value needed for Eq. [12], but makes it easier to see a real flow effect by the non-flatness of 
the distribution in the graph. 

For the correlation between angles determined from two sub-events of different harmonics 
one can write similar relations. For example, the correlation between and ^ is 

(cos(2(*? - **))) = (cos(2(*? - *v)))(cos(2(*£ - * r ))> ( 13 ) 

This expression is proportional to t^i^/er 3 and can be rather small in magnitude, but it can 
be useful for the determination of the relative orientation of the flow of different orders (see 
also the section on Approximations below). On the other hand, the relative orientation also 
can be determined using Eq. (J7|) with k > 1. 

The determination of the event plane resolution. The above relations permit the 
evaluation of the event plane resolution directly from the data. For example, if one knows 
the correlation between two equal multiplicity sub-events (where the resolution of each is 
expected to be the same) then from Eq. fll2]) the resolution of each of them is 

(cos(n(*S. - * P ))> = ^cosW^-v^))), (14) 

where, as before, n = km, and k is not necessarily equal to 1. If the sub-events are correlated, 
then the term inside the square-root is always positive.^ Note that the event plane resolution 



4 For small amounts of flow, fluctuations can cause this term to be negative for even harmonics 



S 



determined in such a way is the event plane resolution of the sub-events. If one wants to use 
the full event (all detected particles from the event) to determine the event plane, then the 
full event plane resolution can be calculated from the sub-event resolution using Eq. (|Tl~D 
or the approximations below, Eqs. (|23| , |24|) , taking into account that the multiplicity of the 
full event is twice as large as the multiplicity of the sub-event. Because Xm = v m /cr is 
proportional to yN, in a case of low resolution where the curves in Fig. [I] are linear, this 
reduces to 

(cos(n(^ m - %))) = >/2<cos(n(tt* - tf r ))). (15) 

If the sub-events are not "equal" , or if one has only correlations between particles in different 
windows, and the resolution in each window can be different, then one needs at least three 
windows to determine the event plane resolution in each of them. In this case, for example, 
the resolution in the first window is determined as |3|,|4| 



<cos(n(<C - tf r ))> 



(cos(n(^ - WJ)) (cos(n(^ - ggj)) 



\ (cos(n(tf*, - 

Approximations. We want to discuss two limits. In a case of weak flow, Xm 1 
(v m <C cr, which also means that the event plane resolution is low) one can expand the 
exponent under the integral in Eq. @ before the integration. It will follow that in this limit 

ex 1 + Jfx m cos(m(tt m -tt P )). (17) 



d(m(^ m - * P )) V 2 

Then the reaction plane resolution can be estimated analytically to be (the case k — 1) 

(cos(m(^ m - %))) w ^ Xm- (18) 

In this case the resolution is linear in x an d, since a is the same for all m, the resolution for 
different m scales with v m . 

Combining Eqs. (|14|)10,18) one can estimate the flow signal directly from the correlation 



between two sub-events 



1W (cos(m(^-^))) 

Vh 2 V N sub ' [U) 



when the total pt is required to be zero. 
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where N SU b is the multiplicity of the sub-events. 

Keeping the second order terms in the expansion of Eq. (Q) it can be shown that 



(cos(2m(^ m - %))) « ^ « - ((cos(m(^ m - * r )))) 2 . 




(20) 



The latter relation for k = 2 is needed as a resolution correction for v 2 in the case where 
the event plane was determined by the directed flow (m = 1). The approximate relations 
(|18|.|(]) are rather accurate for Xm < 0.5. 

Using the above approximations the correlation between flow angles of sub-events of 
different harmonics (m = 1 and m = 2) can be written as 



As stated above, the sign of the left hand side shows the relative orientation of the flow 
of the different harmonics, while the right hand side after the ± sign is always positive. 
Note that all quantities in this equation are evaluated from the data independently, so the 
comparison of the magnitudes is an additional consistency check. 

In the second limiting case of strong flow Xm 3> 1, one can expand the cosine in the 
exponent of Eq. (j^j and get approximately 



which also can be used to calculate the reaction plane resolution. The case of strong flow 
is important mostly for relatively low energy collisions and we will not discuss it in detail. 
Note that in this case (and only in this case) the distribution of m(\l/ m — ty r ) can be described 
by a Gaussian. 



(cos(m(^ m - * r ))) = 0.626657x m - 0.09694x„ + 0.02754^ - 0.002283x™, (23) 
(cos(2m(* m - * r ))) = 0.25 X 2 m - 0.011414x^ - 0.034726;^ + 0.006815xL (24) 



The interpolation formulae are valid for Xm < 3. 

In order to use the above equations to go from the sub-event resolution to the full-event 
resolution, one has to set the equation equal to the sub-event resolution and take the root to 






(22) 



Here we present interpolation formulae of Eq. flTT| ) for the most needed cases of k — 1, 2 
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obtain Xm- One multiplies this by y/2 because Xm is proportional to y/N, and then evaluates 
the full-event resolution from the same equation. Our routines for finding the roots of these 
equations use an iterative method, and the routines for propagating the errors involve a 
calculation for the change in the result with a small finite change in the input. 

Another approximate method is to evaluate Xm — v m / a of the full event from the fraction 



(25) 



of events where the correlation of the planes of the sub-events is greater than tt/2 p6 

N events (m\^ n -^ b m \ > tt/2) _ e-*^/ 4 

Ntotal 2 

This fraction, and therefore the equation, is only accurate when Xm is small. 

III. EVENT-BY-EVENT ANALYSIS OF AZIMUTHAL DISTRIBUTIONS 



Here we study the anisotropies in the azimuthal distributions of particles within a rel- 
atively large rapidity (pseudorapidity) window without determining an event plane. Finite 
multiplicities used in each event for the evaluation of the event flow vector Q n , defined in 
Eqs. causes both finite event plane resolution and also fluctuations in the vector mag- 

nitude Q n . In the case of zero flow (Q n ) = (N). (In this section for simplicity we assume 
W{ = 1.) Anisotropic flow, which shifts the vector Q„ in each event in the (random) flow 
direction by a value v n N, results effectively in broadening the distribution in Q n . Keeping 
the first non-vanishing contribution from flow one has 

(Ql) = (N)+vl(N 2 ). (26) 

This equation can be used to estimate the flow signal v n , which is the average of v n over the 
rapidity region used for the Q n calculation. 

Another method p5| involves fitting the distribution in r n = Q n /N to the theoretical 
formula, which has v n (more exactly its absolute value) as a parameter. In this method the 
distribution of r n has to be fitted by the function 

^exp( __)/ (_ — ), (27) 



r n dr n cr 2 2cr 2 a 2 

where Iq is the modified Bessel function, v n is the parameter of interest, and a is the 
parameter related to finite multiplicity fluctuations. The distribution (p7|) was derived using 
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the central limit theorem, requiring the particle multiplicity iV to be large. Note that in 
principle Eq. fl2"6| ) does not have this limitation, although effectively the "signal to noise 
ratio" is proportional to N and one would need N large to apply the equation to data. 

Note also that both methods do not require the determination of the event plane and 
can be performed within one (pseudo)rapidity window, but they do require this window to 



be rather wide in order to have a relatively large number of particles. The distribution (|27|) 
can be used to select events with larger or smaller values of flow. 

IV. NON-FLOW CORRELATION CONTRIBUTION 

The methods described so far in this paper are correct when correlations induced by flow 
dominate all others. By all others we mean correlations due to momentum conservation J29 



long and short range two- and many-particle correlations (due to quantum statistics, res- 
onances, mini and real jet production, etc.). Below, we discuss the contribution of such 
non-flow correlations to the correlation of the event planes determined from two indepen- 
dent sub-events. Very often the contribution of non-flow correlations scales as 1/N, where 
N is the multiplicity of particles used to determine the event plane. But one should re- 
member that the contribution due to momentum conservation increases with the fraction 
of particles detected, and that the relative contribution of Bose-Einstein correlations can 
be independent of multiplicity (the later can be important if the sub-events are formed by 
random division of all particle from the same phase space region into two subgroups, where 
the particles contributing to the different sub-events can be very close to each other in phase 
space). The study of the effects of non-flow correlations (from the point of view of a flow 
analysis) in real data is rather complicated, so we start with the analysis of Monte-Carlo 
events. 

Non-flow correlations in Monte-Carlo generated events. When the true reaction 
plane is known (as it is in any generated event) the contribution of non-flow correlations 
can be studied by analyzing correlations along the axis perpendicular to the reaction plane 
(y axis). For example, let us consider the correlation between Q a and Q b , the vectors (see 
Eq.([5[)) defined by two independent sub-events. One can think of these vectors as the total 
transverse momentum of all particles of the sub-event (which is the case if the transverse 
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momentum is used as a weight in Eq.fl5|)). Then if there are no other correlations except 
flow 

(Q a Q b ) = (Q a )(Q b ) = (Q a x )(Q b x ), (28) 

the very relation our first method is based on. It was assumed here that the two Q-vectors 
are totally uncorrelated except that both of them are correlated with the reaction plane.0 
If this is not true and there exist other correlations, then their contribution in first order 
would be the same to the correlations of x components and y components. Then 



nan— flow 

{yWy)non-flow = Wx/Wi/ "I" Wj/Vj, 



(Ql) (Ql) + (Q a v Qi)non- f i ow = (Q a x ) (Ql) + (Q a v Qi) (29) 



(Q a x )(Ql)*(QlQl)-{Q a y Q b y)- (30) 

Analyzing the correlation {QyQ b y ) (and in particular (sin(n(\I/" — \I/ r )) sin(n(\l/^ — ^ r ))) one 
can estimate the corrections to the formulae (8,12-15). 

Non-flow correlations in real data. The direct application of the above described 
method to real data is not possible. What can be done is the analysis of similar correlations 
using, instead of \l/ r , the event plane derived from the second harmonic, where as the analysis 
of different models shows, the contribution of non-flow effects is significantly less. Then with 
the (second harmonic) event plane resolution known one can carry out the above analysis. 

There exists in the literature other methods for estimating and accounting for non- 
flow correlations, see, for example refs. |]30|j33| . Here we briefly describe the method []33 



which was applied to the data of the WA93 Collaboration [[34]]. It was proposed to 
characterize the non-flow correlation contribution by the value of the parameter 

c _ (Q a Q b ) - (Q a x )(Qi) „ (Q a Q b ) - mm m 

v / ((Q a ) 2 ) v / ((Q fe ) 2 ) ~ N 

where iV is the sub-event multiplicity and, as in section |T|, for simplicity we assume Wi — 1. 
The parameter c can strongly depend on the particular choice of sub-events, but if the 



5 Strictly speaking there is one more assumption here, namely that the strength of the flow does 
not fluctuate event by event. 
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non-flow contribution is dominated by two-particle correlations, it is largely independent of 



multiplicity p3||. The non-flow correlation changes the distribution of event flow vectors, 



and in particular Eq. (|26|) now reads 

(Ql)K(N)+v 2 n (N 2 ) + 2c(N), (32) 

If the parameter c is relatively large, it can bias the results derived from an application 
of Eq. fl2"ED to the data. However, the non-flow correlations contribute to the shape of the 



distribution in Q n (Eq. (^)) in a different way than flow. They produce mostly a change 
in a, the "Gaussian width" of the distribution (Eq.(|T0"D), which is modified to 

a 2 = ^(l + 2c), (33) 

and the parameter c can be directly extracted from the data. The analysis |34| of S+S and 



S+Au data gave (for their particular sub-event selection) c ~ 0.034 ± 0.025 which is not 
negligible compared with t> 2 ~ 0.04-0.05 found in this analysis at multiplicities of about 
30-50. Another possible consequence of non-flow contributions could be the change in the 
shape of the distribution of the difference of flow angles of sub-events [[331 . 



V. SIMPLE WAY TO INTRODUCE FLOW IN A MONTE-CARLO EVENT 

GENERATOR 

Sometimes in order to investigate different detector effects or the reliability of the method, 
one needs to introduce anisotropic flow into a Monte-Carlo event generator. It can be done 
by changing the azimuthal angle of each particle (and consequently changing the density in 
the azimuthal angle space) 

4> — 4>' = 4> + A<f>, (34) 

where 

_2 

A(f) = Y, — v n sm(n((/)-ip )), (35) 

n n 

v n are the n parameters of the transformation, and ipo is the direction of the added flow. 
The v n can be functions of rapidity and transverse momentum and, in particular, the v n for 
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n odd should reverse sign in the backward hemisphere. One can check that such a change 
in the azimuthal angle results in the required change in the distribution. To first order in 
v n , v n = (cos(n(0' — \? r ))) = v n . Small higher order corrections (~ v*), if needed, can also 
be easily calculated, for example, numerically. 

VI. DISCUSSION 

Higher harmonics (n > 3). Note that the flow analysis methods presented in this 
paper are valid for all harmonic orders of anisotropy. The higher harmonics look at the event 
with higher symmetry on a finer scale. It should be emphasized that the study of anisotropic 
flow corresponding to harmonic orders of n > 3 has interesting aspects. For instance, one 
would expect large differences between the theoretical predictions of hydro and cascade 
type models in the higher harmonics of the particle azimuthal distributions. Also there 
are physics processes which could lead to non-zero higher harmonics. It is widely accepted 
that one of the main reasons for pion "anti-flow" , directed flow in the direction opposite to 
that of the protons, is caused by pion shadowing by co-moving nucleons. Such shadowing 
could affect pion higher harmonic azimuthal distributions, and affect them differently at 
different rapidities. Another effect is the transition from out-of-plane elliptic flow (squeeze- 
out) which is very important at low energies to in-plane elliptic flow which is the main effect 
at high energies. At the transition beam energy both effects may be important and the 
fourth harmonic may peak at this energy. Even at high beam energies, the out-of-plane 
squeeze-out effect could dominate at short times and the in-plane expansion at long times, 
leading to an overall fourth order harmonic coefficient. 

Transverse radial and anisotropic flow. We would like also to emphasize that the 
above methods permit the reconstruction of the triple differential distribution with respect 
to the reaction plane, and in particular are convenient for the analysis of the pt dependence 
of the anisotropy. The importance of such a study was stressed in ||35|| , where the effect 
of the interplay of (transverse) radial and directed flow has been studied. In that analysis 
transverse directed flow is considered a result of the movement of an effective source in 
the transverse plane. It is assumed that in the source rest frame the first moment of the 
azimuthal distribution (vi) is zero, and the final anisotropy appears only as a consequence 
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of the source movement in the transverse direction. 

In this case the transverse momentum dependence of v\ has a rather specific shape. The 
radial expansion results in decreasing v% at low p t . For some sets of parameters (temperature, 
radial and directed velocities) it becomes negative. Physically it corresponds to the case, 
when particle production with such a value of pt is more probable from the part of the 
effective source which moves in the opposite direction from the flow direction, because in 
this case the directed flow and the radial flow tend to compensate each other. 

The same considerations can be applied to elliptic flow as well. If one assumes that 
elliptic flow is a consequence of more rapid expansion of the effective source in some plane, 
then the p t dependence of v 2 would exhibit exactly the same features as has been observed 
for Vi in the case of directed flow. 

Pair-wise azimuthal correlations. Two-particle large-angle azimuthal correlations 



are often proposed as a tool to study anisotropic flow [36]. Not rejecting this possibility we 
note that the expected signal in such correlations can be very small in magnitude. It can be 
easily shown that the pair-wise distribution in the azimuthal angle difference (A0 = 0i — ^2) 
is 

dN pairs 



oc (1 + ]T 2vl cos(nA0)). (36) 



Note that the "signal" is v^, and thus is small. It does not imply though, that the "signal 
to noise ratio" is small, and the method, in principle, can be successfully applied to data 
to obtain the flow signal. However the reconstruction of the triple differential distribution 
with respect to the reaction plane (the goal of the flow analysis) becomes more involved. 

This method does not require the determination of the event plane. Usually such a 
distribution is constructed using all possible pair combinations in each event. Note that 
in this case each particle enters into many pairs (on the order of the mean multiplicity 
of an event) and, consequently, the pairs are not statistically independent. Thus special 
precautions have to be taken for the evaluation of the error of the results [pT| . 
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FIG. 1. The event plane resolution for the n— th (n = km) harmonic of the particle distribution 
with respect to the m— th harmonic plane, as a function of Xm = v m /a. 



20 



